Prioritizr Runs

We used the fPrioritizrRun() function to streamline the runs. fPrioritizrRun() arranges all the data and makes sure that they’re in the right format, then runs problem() and solve(). This function also saves the summaries each plan (cost, % area protected and number of representative features that achieved 30% or more protection).

The relevant scripts for this step is found in ~/scripts/ with the prefix 10.

# the function is in:
source("scrpits/10a_Prioritizr_fPrioritizr.R")
# to run the function with climate-smart data
source("scripts/10b_Prioritizr_fPrioritizrAQMSmartRun.R")
# to run the function with climate-uninformed data
source("scripts/10d_Prioritizr_fPrioritizrAQMUninformedRun.R")
# to create no-regret closures with the climate-smart plans

To create no-regret plans, we used thre function fCreateNoRegret() which intersected the selected PUs in all three climate scenarios. The relevant scripts for this step is found in ~/scripts with the prefix 11.

# the function is in:
source("scripts/11a_NoRegret_fCreateNoRegret.R")
# the runs are found in:
source("scripts/11b_NoRegret_fCreateNoRegretAQMRuns.R")

As mentioned in the previous section, we will only be showing results of Runs 2 (Max = 0.9, Min = 0.1), 5 (Max = 0.6, Min = 0.1) and 8 (Max = 0.3, Min = 0.1).

Run 2: Maximum Target = 90%

Climate-smart runs

SSP1-2.6 run

SSP126_target90 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/02_Target90/solution_SSP126.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario    cost perc_area_protected protected_PU
##    <dbl> <chr>      <dbl>               <dbl>        <dbl>
## 1     90 SSP1.26  267943.                26.0           59

SSP2-4.5 run

SSP245_target90 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/02_Target90/solution_SSP245.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario    cost perc_area_protected protected_PU
##    <dbl> <chr>      <dbl>               <dbl>        <dbl>
## 1     90 SSP2.45  267309.                26.9           61

SSP5-8.5 run

SSP585_target90 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/02_Target90/solution_SSP585.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario    cost perc_area_protected protected_PU
##    <dbl> <chr>      <dbl>               <dbl>        <dbl>
## 1     90 SSP5.85  313651.                29.8           60

No-regret plan

noregret_target90 <- readRDS("outputs/11_NoRegret/11b_AQMRuns/noregretclosures_Target90.rds")

Climate-uninformed run

uninformed_target90 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/02_Target90/solution_uninformed.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario      cost perc_area_protected protected_PU
##    <dbl> <chr>        <dbl>               <dbl>        <dbl>
## 1     90 uninformed 744951.                73.5           66

Run 5: Maximum Target = 60%

Climate-smart runs

SSP1-2.6 run

SSP126_target60 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/05_Target60/solution_SSP126.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario    cost perc_area_protected protected_PU
##    <dbl> <chr>      <dbl>               <dbl>        <dbl>
## 1     60 SSP1.26  132371.                17.3           54

SSP2-4.5 run

SSP245_target60 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/05_Target60/solution_SSP245.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario    cost perc_area_protected protected_PU
##    <dbl> <chr>      <dbl>               <dbl>        <dbl>
## 1     60 SSP2.45  131116.                18.0           54

SSP5-8.5 run

SSP585_target60 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/05_Target60/solution_SSP585.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario    cost perc_area_protected protected_PU
##    <dbl> <chr>      <dbl>               <dbl>        <dbl>
## 1     60 SSP5.85  155441.                19.9           55

No-regret plan

noregret_target60 <- readRDS("outputs/11_NoRegret/11b_AQMRuns/noregretclosures_Target60.rds")

Climate-uninformed run

uninformed_target60 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/05_Target60/solution_uninformed.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario      cost perc_area_protected protected_PU
##    <dbl> <chr>        <dbl>               <dbl>        <dbl>
## 1     60 uninformed 369968.                49.0           63

Run 8: Maximum Target = 30%

Climate-smart runs

SSP1-2.6 run

SSP126_target30 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/08_Target30/solution_SSP126.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario   cost perc_area_protected protected_PU
##    <dbl> <chr>     <dbl>               <dbl>        <dbl>
## 1     30 SSP1.26  49676.                8.69           39

SSP2-4.5 run

SSP245_target30 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/08_Target30/solution_SSP245.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario   cost perc_area_protected protected_PU
##    <dbl> <chr>     <dbl>               <dbl>        <dbl>
## 1     30 SSP2.45  49934.                9.01           39

SSP5-8.5 run

SSP585_target30 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/08_Target30/solution_SSP585.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario   cost perc_area_protected protected_PU
##    <dbl> <chr>     <dbl>               <dbl>        <dbl>
## 1     30 SSP5.85  57173.                9.96           38

No-regret plan

noregret_target60 <- readRDS("outputs/11_NoRegret/11b_AQMRuns/noregretclosures_Target60.rds")

Climate-uninformed run

uninformed_target30 <- readRDS("outputs/10_Prioritizr/10b-c_AQMRuns/08_Target30/solution_uninformed.rds")

cost is in USD
perc_area_protected is the percentage of the Pacific ABNJ protected
protected_PU is the number of features that have >= 30% (by area) protection

## # A tibble: 1 x 5
##   target scenario      cost perc_area_protected protected_PU
##    <dbl> <chr>        <dbl>               <dbl>        <dbl>
## 1     30 uninformed 129125.                24.5           42